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Counterion distributions and effective electrostatic interactions between spherical macroions in 
polyelectrolyte solutions are calculated via second-order perturbation (linear response) theory. By 
modelling the macroions as continuous charge distributions that are permeable to counterions, an- 
alytical expressions are obtained for counterion profiles and effective pair interactions in solutions 
of star-branched and microgel macroions. The counterions are found to penetrate stars more eas- 
ily than microgels, with important implications for screening of bare macroion interactions. The 
effective pair interactions are Yukawa in form for separated macroions, but are softly repulsive and 
bounded for overlapping macroions. A one-body volume energy, which depends on the average 
macroion concentration, emerges naturally in the theory and contributes to the total free energy. 



I. INTRODUCTION 

Polyelectrolytes (PEs) are ionizable polymers that 
dissolve in a polar solvent, such as water, by 
dissociating into polyvalent macroions and small 
oppositely-charged counterions 0. Electrostatic 
interactions between macroions, mediated by sur- 
rounding microions (counterions and salt ions), con- 
tribute to the unique macroscopic properties of PE 
solutions, which are the basis of many industrial 
applications involving polymer- water systems 0|- 
Common synthetic examples of PEs are polyacrylic 
acid, used in gels and rheology modifiers, and 
polystyrene sulfonate, a component of reverse osmo- 
sis membranes. Naturally occurring examples are 
biopolymers, such as DNA, proteins, and starches. 
Colloidal in size, PEs are also routinely added as 
flocculants and stabilizers to colloidal suspensions, 
such as foods and water-based paints % 11. De- 
pending on PE concentration, adsorption or grafting 
of PE chains onto surfaces of colloidal particles can 
either induce flocculation, by bridging particles, or 
impart electrosteric stabilization. 

Conformations of PE macroions and electrostatic 
interactions between macroions are strongly influ- 
enced by the distribution of microions. If dis- 
persed in solution, microions act to screen the bare 
Coulomb interactions between ionized monomers. If 
condensed on the macroion chains, microions may 
reduce the macroion charge j^. Linear PE chains 
whose monomers are sufficiently weakly interacting 
- either because weakly charged or because of strong 
microion screening or condensation ~ may form 
random-walk coils with roughly spherical conforma- 
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tion. With increasing charge and screening length, 
linear chains stretch into non-spherical conforma- 
tions because of electrostatic repulsion between ion- 
ized monomers Q, • The extent of elongation de- 
pends on the chain charge density, salt concentra- 
tion, and solvent quality. Highly charged chains 
in good solvents {e.g., DNA in water) often form 
stiff rod-like macroions, whose effective interactions 
and complex phase behavior (such as bundling) have 
been widely studied [1, I3- In poor solvents, suf- 
ficiently highly charged chains may form necklaces 
of compact globules joined by narrow threads, as 
predicted by theory [10|| and confirmed by simula- 
tion [ni. 

Although many common PEs are linear, other 
topologies can be readily synthesized. Examples 
are stars, microgels, micelles, and brushes. Star 
polymers consist of chains chemically grafted 
or adsorbed to a common microscopic core. Micro- 
gels are mesoscopic polymer networks^ynthesized 
by polymerization in microemulsion Micelles 
are formed by association of charged diblock (am- 
phiphilic) copolymers Brushes are forme d by 

grafting PE chains onto a mesoscopic solid core 'T^ . 
Solutions of spherical stars, microgels, micelles, and 
brushes can be regarded as colloidal suspensions of 
soft macroions that are permeable to microions. 

Electrostatic interactions in charged colloids have 
received much attention in recent years lUl, mo- 
tivated largely by anomalous phase behavior that 
is unexplained by the classic Deriaguin-Landau- 
Verwey-Overbeek (DLVO) theory [T3- Most studies 
have been restricted, however, to hard, impermeable 
macroions. The objectives of this paper are first, 
to explore implications of microion penetration for 
screening of effective electrostatic interactions be- 
tween spherical macroions, and second, to lay a foun- 
dation for future studies of thermodynamic phase 
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behavior of PE solutions. Our approach is based on 
a recently-proposed theory of effective interactions 
in charged colloids 18], which we adapt here from 
hard to penetrable macroions and apply to spherical 
star-branched and microgel macroions. 

The remainder of the paper is organized as fol- 
lows. Section ^ describes the assumed model of 
PE solutions. Section ITTTl reviews the theoretical ap- 
proach, based on second-order perturbation (linear 
response) theory. Sections IIVI and present ana- 
lytical and numerical results for counterion profiles 
and effective interactions in bulk solutions of star 
and microgel macroions. Finally, Sec. I VII closes with 
a summary and conclusions. 



II. MODEL 

Adapting the primitive model of ionic liq- 
uids ilS], the model system comprises Nm spheri- 
cal macroions, of radius a (diameter a = 2a) and 
charge —Ze, and Nc point counterions of charge ze 
dispersed in an electrolyte solvent in a volume V 
at temperature T. Assuming, for simplicity, a sym- 
metric electrolyte and equal salt and counterion va- 
lences, the electrolyte contains Ng point salt ions 
of charge ze and Ns of charge —ze. The microions 
thus number N+ = Nc + Ng positive and N- = Ng 
negative, for a total of N/^^ = Nc + 2Ns. Global 
charge neutrality in a bulk solution constrains av- 
erage macroion and counterion number densities, 
n-m — Njn/V and Uc — Nc/V, via Zrim — zuc- The 
polar solvent is treated as a continuum, character- 
ized by dielectric constant e that acts only to reduce 
Coulomb interactions between ions. 

The local number density profiles of counterions, 
Pcir), and of macroion monomers, Pmon(?'), are mod- 
elled as spherically symmetric, continuous distri- 
butions. Spherical symmetry is a reasonable ap- 
proximation, considering that equilibrium averaging 
over macroion orientations tends to smear out any 
anisotropy. Furthermore, discreteness of the charge 
distributions can be ignored if we restrict considera- 
tion to length scales exceeding the scale of discrete- 
ness. 

In general, the counterions are distributed over 
three regions: (1) the immediate vicinity of the PE 
chains making up the macroions, (2) the region in- 
side of the macroions but away from the chains, and 
(3) the region outside of the macroions. Counterions 
in the first two regions are trapped by the macroions, 
while those in the third region are free. Within the 
first region, the counterions may be either condensed 
on a chain or free to move along a tube surrounding 
a chain. These chain- localized counterions, whether 



condensed or mobile, tend to distribute uniformly 
along the chains to favor local charge neutrality. In 
our model, counterions in region (1) simply renor- 
malize the effective macroion valence Z . 

The detailed form of the monomer density pro- 
file depends on the macroion conformation. For 
star-branched macroions. Coulomb repulsion be- 
tween charged monomers tends to stiffen and radi- 
ally stretch the chains into a porcupine conforma- 
tion [23 . We assume the ideal case of fully stretched 
chains and model the monomer density profile by 
PmoniT 1^ a) — Z / {Anar"^) , where r is the radial dis- 
tance from the star's center. For microgel macroions, 
the dense network of chains is well approximated by 
a uniform monomer distribution, and is modelled 
here by Paion{r < a) — 'SZ/ {Ana'^). This distribu- 
tion may also approximate a weakly-charged linear 
PE chain with a spherical random-coil conformation, 
although a Gaussian distribution may then be more 
accurate. For both the star and microgel models, 
the monomer density profile is cut off sharply at the 
macroion surface: Pmonir > a) = 0. 



III. THEORY 

For the model PE solutions described above, the 
theoretical challenge is to predict the distributions 
of microions inside and outside of the macroions and 
the effective interactions between macroions. Fol- 
lowing the same general strategy as applied pre- 
viously to charged colloids 0, |^, we reduce 
the multi-component mixture to an equivalent one- 
component system governed by effective interac- 
tions, which are approximated via perturbation the- 
ory. For clarity of presentation, we initially ignore 
salt ions. The Hamiltonian then decomposes into 
three terms: 

H = H„,{{R}) + Hc{{r}) + H„,c{m, {r}), (1) 

where {R} and {r} denote collective coordinates of 
macroion centers and counterions, respectively. The 
first term, 

is the bare Hamiltonian for macroions with kinetic 
energy Km that interact via the bare pair potential 
Vrnmi^) ^t centcr-centcr separation r. The form of 
Vmm{T) depends on the macroion conformation and 
is specified in the Appendix. The second term in 
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Eq. O, 

Hc = K, + - ^ «,e(|r,-r,|), (3) 

is the Hamiltonian for counterions with kinetic en- 
ergy Kc that interact via the Coulomb pair potential 
Vcc{t) = z'^e^/er. The third term in Eq. (fT|l. 

1=1 j=i 

is the macroion-counterion interaction. For spherical 
macroions, 

vmcir) = I (5) 
[ w<(r), r < a, 

where the interaction inside a macroion, i'<(r), de- 
pends on the macroion conformation and is specified 
in Sec. IIVI For later reference, we note that Eq. Q 
also may be expressed in the form 



Hmc = J dR pm{R) J drpc{r)v„ 



|R-r|), (6) 



where pm (R.) 



Rj) and Pc(r) 



the macroion and counterion mmr- 
ber density operators, respectively. 

The mixture of macroions and counterions is for- 
mally reduced to an equivalent one-component sys- 
tem by tracing over counterion coordinates. Denot- 
ing counterion and macroion (classical) traces by ( )^ 
and ( respectively, the canonical partition func- 
tion can be expressed as 

Z = ((exp(-/3i7)),)„^ = (cxp(-/^i/eff)),„ , (7) 

where Hcs = H„i+Fc is the effective one-component 
Hamiltonian, (3 = l/ksT, and 



(8) 



F, = -fcsTlnfexp -/3(i/c + i?mc) 



is the free energy of a nonuniform gas of counterions 
in the presence of the macroions. 

At this stage, approximations are necessary for 
the counterion free energy. It is first convenient to 
convert the counterion Hamiltonian to the Hamilto- 
nian of a classical one-component plasma (OCP) of 
counterions by adding to He, and subtracting from 
iJmcj the energy of a uniform compensating nega- 
tive background Eb = —NcncVcc{0)/'2, where 
Vcc{0) is the /c — > limit of the Fourier transform of 



Vcc{r). Now regarding the macroions as an "exter- 
nal" potential for the OCP, we invoke perturbation 
theory [Hill 113 and write 



Focp 



{H',nc)\ ■ 



(9) 



where Fqcp = -fcsTln (cxp[-/3(i/, + ^b)])^ is 
the OCP free energy, the A-integral charges the 
macroions, -ff^ic = ^mc — represents the perturb- 
ing potential of the macroions acting on the counte- 
rions, and (H^e)^ is the mean value of this potential 
in a solution of macroions charged to a fraction A of 
their full charge. Further progress is facilitated by 
expressing Hmc [Eq. lO] in terms of Fourier compo- 
nents: 

+ ■^lim[t)™c(fc)/5,„(-k)(p,(k))J, (10) 

V k—>0 

where Vmc{k) is the Fourier transform of Eq. (jsj 
and where /5m(k) = J2j=i exp(ik • R^) and /5c(k) = 

^^-^ exp(ik • Tj) are Fourier components of the 
macroion and counterion densities. 

In first-order perturbation theory, the response of 
the counterion plasma to the macroions is ignored. 
Here we apply second-order perturbation (linear re- 
sponse) theory, in which the counterions are as- 
sumed to respond linearly to the macroion external 
potential: 



Pc(k) = x(fc)^'mc(fc)/Om(k) 



fc^O, (11) 



where x(fc) is the linear response function of the 
OCP. Note that the k ^ limit here, and in 
Eq. l(Tn|) , must be treated separately, since the aver- 
age counterion density, ric = /5c(0), does not respond 
to the macroion charge, but rather is fixed by the 
constraint of global charge neutrality. 

Upon combining Eqs. the effective 

Hamiltonian can be recast in the form of the Hamil- 
tonian of a pairwise-interacting system: 



^cff — Kn 



IE.. 



Kc + t: y . I'cffdR^-R,!) + Eo, 



(12) 

where Vesir) (r) + Vind{r) is an effective 

macroion pair interaction that combines the bare 
macroion interaction with a microion-induced inter- 
action 



Wind(fc) = X{k) [Vmc{k)Y 



(13) 
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The final term in Eq. H12|l is the volume energy, for- 
mally given by 



En — F( 



OCP 



N„i lim 



— Imi Wi„d(?-) 

-^Wind(fc) + ricVmcik) + ^^Vcc{k) 



(14) 



which is a natural by-product of the one-component 
reduction. Although independent of the macroion 
coordinates, the volume energy depends on the av- 
erage macroion density and thus can influence ther- 
modynamics. 

The linear response function is proportional to the 
corresponding static structure factor, S(k), which 
may be obtained from liquid-state theory In 
practice, the OCP is weakly correlated, with cou- 
pling parameter F = Xb/olc ^ 1, where \b — 
I3e'^/e is the Bjerrum length and Oc = (3/47rnc)^/'^ 
is the counterion sphere radius. For example, for 
macroions of diameter a = 100 nm, valence Z — 100, 
and volume fraction ry = (7r/6)ri,„(T'^ = 0.01, in 
water at room temperature {\b = 0.714 nm) , w e 
find F ~ 0.014. As for charged colloids [11 1^, 
we adopt the random phase approximation (RPA), 
which is accurate for weakly-coupled plasmas. The 
RPA equates the two-particle direct correlation 
function of the OCP to its exact asymptotic limit: 
c(2)(r) = —Pvccif). Using the Ornstein-Zernike re- 
lation, S{k) = 1/[1 — ncCp'\k)\, the linear response 
function then takes the analytical form 



X(fc) = ~l3n,S{k) = - 



(l + K2/fc2)^ 



(15) 



where k — v^tto^Jz^Ab is the inverse Debye screen- 
ing length. Note that since permeable macroions 
do not exclude counterions from their interiors, the 
excluded-volume corrections required for hard col- 
loidal macroions |3 ^-re not relevant here. With 
x{k) specified, the counterion density can be ex- 
plicitly determined from Eqs. © and ifTT)) for a 
given macroion distribution (see Sec. IIV|I . Finally, 
salt is easily introduced via additional microion re- 
sponse functions. In the process, the pair inter- 
action and volume energy are unchanged, except 
for a redefinition of the screening constant as k = 
^47r(nc + 2ns)z'^\B, where Us is the average num- 
ber density of salt ion pairs. 

It is worth noting the formal equivalence of 
the present theory to linearized Poisson-Boltzmann 
(DLVO) theory. Both are mean-field theories in the 
sense that they ignore fluctuations in microion dis- 
tributions. An advantage of linear response theory, 
however, is that it encompasses the volume energy. 



which can be important for describing phase behav- 
ior milillllililiil. Moreover, response theory 
can be straightforwardly generalized to incorporate 
nonlinear response, which entails both many-body 
effective interactions and corrections to the pair po- 
tential and volume energy p^ . In contrast, non- 
linear Poisson-Boltzmann theory is practical only 
for the simple boundary conditions afforded by cell 
models. For simplicity, higher-order nonlinear ef- 
fects are here ignored. 

Equations Hll(l - (|14|1 constitute the main for- 
mal expressions from linear response theory. Ex- 
plicit calculations require specifying the counterion- 
macroion interaction u< (r) in Eq. Q for specific 
macroion models. Below, we apply the theory to ob- 
tain analytical and numerical results for counterion 
profiles and effective interactions in bulk solutions 
of spherical star-branched and microgel macroions. 



IV. ANALYTICAL RESULTS 



Star Macroions 



For our idealized model of a star-branched 
macroion with l/r^ monomer density profile. 
Gauss's law gives the electric field as 

r Ze 



E{r) = 



ear ' 



r > a 
r < a. 



(16) 



Integration over r yields the electrostatic potential 
energy between a star and a counterion: 



Vmc{r) = 



Zze^ 



r > a 



[l-ln(£)], r<a, 



whose Fourier transform is 



(17) 

sinc(fca), (18) 



with sinc(a;) = Jq du sm{u)/u. We can now calcu- 
late the counterion number density around a single 
macroion in the dilute limit, where /0„i(k) = 1. From 
Eqs. 1)11(1 . I|15|) . and (|18|l . the Fourier component of 
the counterion density profile is given by 



Pcik) = 



Z 



z ka{k^ + K^) 
whose real-space form is 
Z K 



sinc(/ca), (19) 



Pc{r) 



z Airar 



sinhc(Ka) e 



r > a 



(20) 
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Z K 

Pc(r) = — [Ec(Ka, Kr) + 2 sinhc(/i 



z Sirar 

Ec(-Ka, -Kr) e''^ r < a, (21) 



where 



sinhc(a;) = 



du 



sinh(u) 



(22) 



and 

Ec(a;i,a;2) = / 



du ■ 



In 



(23) 

which can be efhciently computed from the first few 
terms of the rapidly converging series expansions. 
Approaching the macroion center, the counterion 
density profile varies more gradually than the 
macroion monomer density profile, diverging loga- 
rithmically, according to 



lim pc{r) 



Ana 



1 - In 



(D 



(24) 



Integrating Eq. 1)21(1 over the spherical volume of the 
macroion yields the fraction of countcrions inside a 
star: 



fir 



Z . 

—4tt 
Z 



1 



dr r^pc{r) 



1 + — ) e""" sinhc(Ka). (25) 



Note the clear predictions that (1) the counterion 
distribution is determined entirely by ko, or the di- 
mensionless ratio of the macroion radius and the De- 
bye screening length, and (2) the fraction of countc- 
rions inside increases monotonically with na. Thus, 
for fixed macroion radius, /in increases with increas- 
ing macroion valence and concentration. This re- 
sult is physically sensible: the shorter the screening 
length, the shorter the range of the counterion re- 
sponse, and thus the tighter the localization of coun- 
tcrions around the macroion centers. 

From Eqs. ((T^ and ((TH|l . the induced electrostatic 
pair interaction is given by 



t'ind(fc) 



Fourier transforming, we obtain 



sinc^(A;a). 



(26) 



mn^Z^e^K^a^ 



dx ■ 



sin{xr/a) 



X smc X. 



(27) 



For nonoverlapping stars, Eq. H27() can be reduced 
to the analytical form 



Vind{r > 2a) 



Z'e' Z'e' 



sinhc(Ka) 



Ka 



r 

(28) 



Since nonoverlapping macroions interact via a bare 
Coulomb potential, «,„,„(?■) — Z^e^ /er, the effective 
pair interaction for this case is 



2„2 r 



Vcs{r > 2a) = 



Z'e 



sinhc(Ka) 



(29) 



Thus, at the level of linear response, nonoverlap- 
ping star macroions interact via an effective Yukawa 
(screened-Coulomb) pair potential. The screening 
constant, k, in the potential depends on the to- 
tal density of microions - inside and outside of 
the macroions - since all microions respond to the 
macroion charge. Note that the potential has the 
same r-dependence as the DLVO potential for hard 
colloidal macroions 0, 0| , 



vuiMo{r) 



exp(Ka) 
I + Ka 



, r > 2a, 

r 

(30) 

differing only in the macroion-size-dependent ampli- 
tude. For overlapping stars, the bare macroion in- 
teraction is somewhat more complex and is relegated 
to the Appendix. 

Finally, from Eqs. (HH), (HHJl, (EHl), and (|13), the 
volume energy is obtained as 



sinc^a; 



e Jq x'^yx'^ + K^a^) 

(31) 

For weakly-coupled microion plasmas, the OCP free 
energy may be approximated by its ideal-gas limit: 

i^ocp = TV+H^+A^) - 1] + iV_[ln(n_A3)-l], 

(32) 

where A is the thermal wavelength. Note that if Z 
is allowed to vary (with counterion condensation), 
then the volume energy per macroion must be aug- 
mented by the self (Hartree) energy of a macroion: 
Uh = Z'^e^jta. The first term in Eq. (EJ^ long 
recognized as important for phase behavior [J , rep- 
resents the entropy of free counterions; the sec- 
ond term accounts for the cohesive electrostatic en- 
ergy of microion-macroion interactions. The vol- 
ume ene rgy, analogous to its counterpart for charged 
colloids [isL I23. I25I |25| . depends on the average 
macroion concentration and thus has the potential to 
influence phase behavior and other thermodynamic 
properties. Equations (HOJl, (EU, (EHl, (EHl, and 
are the main analytical results for star macroions. 
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It is important to emphasize that the present ap- 
proach, while including the entropy of the counte- 
rions, neglects the configurational entropy of the 
macroions by assuming rigid (fully stretched) PE 
chains. Recently, Jusufi et al. f28\ modelled pair in- 
teractions between PE stars by both molecular dy- 
namics simulation and a variational free energy that 
incorporates chain flexibility. An important con- 
clusion of their study is that pair interactions are 
dominated by counterion entropy. Our approach 
is complementary: while the macroion model ne- 
glects chain flexibility, which is reasonable at least 
for nonover lapping stars, the linear response theory 
refines somewhat the modelling of the counterion 
distribution. 



B. Microgel Macroions 

For our model of microgel macroions, we apply 
exactly the same procedure as in Sec. IIV Al The 
electric field of a uniformly-charged sphere is 



E{r) 



Ze 



r > a 
r < a, 



(33) 



which integrates to give the macroion-counterion in- 
teraction, 



Vmcir) 



er ' 



r > a 
r < a. 



(34) 



Equation H34|) Fourier transforms to 



Vl-nZze^ \ sin(fca) 
cos(fca) 



ka 



(35) 



which, when substituted into Eq. yields the 

Fourier transform of the counterion density profile 
around a single macroion. 



/5c(fc) = - 



z k'^a'^ik' + k'^) 



cos(fca) 



sin(fca) 



(36) 

The Fourier transform of Eq. H3()|l gives the real- 
space counterion density profile. 



Pc{r) 



Pc{r) = 



Z 3 
z Aira^r 



Z 



sinh(Ka) 
cosh(Ka) a 



r > a 



(37) 



- 1 



a 



r < a, 



' sinh(Kr) 



(38) 



which approaches a constant as r 
Z 3 
z Aira^ 

Integration of Eq. H38() yields an analytical result 
for the internal fraction of counterions: 



Pcir - 0) 



[l - (1-1- Ka)e'' 



(39) 



/i. 



1- 



3 

Ka 



1 



1 

Ka 



cosh(Ka) 



sinh(Ka) 
Ka 

(40) 

Again the theory predicts a counterion distribution 
depending only on the ratio of macroion radius to 
screening length and an internal counterion fraction 
that increases monotonically with this ratio. This 
prediction may be compared with that of Oosawa's 
"two-phase" approximation T|, which assumes uni- 
form (but differing) counterion concentrations inside 
and outside of the macroions. According to the lat- 
ter approach, for spherical macroions with volume 
fraction 77, the condition for equilibrium between free 
and bound counterions, in the absence of salt ions. 



In 



= In 



V 



, f , -.^ . rZ^{l-h^){l-v"% 

(41) 

which must be solved numerically for /in. The pre- 
dictions of Eqs. (|40|l and H41I) arc compared below 
in Sec. El 

Next, substituting Eq. into Eq. the in- 

duced pair interaction is 



Vind(fc) 



367rZ2e2 



cos(fca) 



/fc6a4(fc2 + K^) 
sm{ka) 



ka 



(42) 



For nonoverlapping macroions, the Fourier trans- 
form of Eq. (|42|l is straightforward to evaluate and 
yields an effective pair interaction 



Z'e' 



9 



cosh(Ka) 



sinh(Ka) 



Ka 



r > 2a. 



(43) 



As for star macroions, a Yukawa form is predicted, 
but with a different amplitude. The case of overlap- 
ping macroions is left to the Appendix. 

Finally, from Eqs. (dj), and the volume 
energy is obtained as 

3Z2e2 



Eo 



Focp 
3 

AK^a^ 
{N+- 



N„ 



1 
2 

Ka K 



N-) 



ksTK^a' 



(44) 
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Equations gOl), and are the 

main analytical results for microgel macroions. 



V. NUMERICAL RESULTS 

The theory developed above can be applied to 
solutions of arbitrary ionic strength, under the 
assumption that the macroion PE chains remain 
stretched. In order to highlight the role of the coun- 
terions, we present numerical results for salt-free so- 
lutions. Within the model considered, the effect of 
salt is merely to increase the Debye screening con- 
stant. Furthermore, we consider the case of mono- 
valent counterions (z = 1) in aqueous solutions at 
room temperature {Xb = 0.714 nm). Figure ^ il- 
lustrates the form of the counterion number den- 
sity profiles [Eqs. (HJ, and (j2Z|l, inside 
and outside of a macroion. Inside a star macroion 
the counterion density diverges logarithmically to- 
wards the center, while inside a microgel macroion 
Pcir) remains finite. Evidently, counterions pene- 
trate stars more easily than they do microgels. This 
property is also reflected in the internal counterion 
fractions [Eqs. H25|l and (|40|) ]. functions of na only, 
which are shown in Fig. In Fig. we com- 
pare predictions of linear response theory [Eq. 
with those of Oosawa's two-phase approximation 
[Eq. H41|l ] for uniformly-charged spherical (microgel) 
macroions 0] . Both approaches qualitatively predict 
an increase in the fraction of bound counterions with 
increasing macroion concentration. However, linear 
response theory predicts a considerably more grad- 
ual accumulation of bound counterions than does the 
two-phase approximation. 

Counterion penetration strongly influences screen- 
ing of bare macroion interactions. The effective pair 
potentials, Vcs{r), and corresponding forces, F{r) = 
— dvefi(r)/dr, are shown in Figs. |31 and 0] respec- 
tively. Beyond overlap the effective interaction has 
Yukawa form, with the amplitude depending on the 
type of macroion. Figure[Slcompares the variation of 
the macroion-size-dependent amplitude of v(r > 2a) 
with the Debye screening constant for the two per- 
meable macroions and for hard macroions. Evi- 
dently, the greater the permeability of the macroions 
to counterions, the weaker the amplitude of long- 
range repulsion. For overlapping macroions, the 
bare charge distribution combined with counterion 
penetration leads to softly repulsive interactions. 
Note that the interactions are bounded: they do not 
diverge as the macroions approach complete overlap. 
It must be emphasized that the effective interactions 
presented in Figs. 13 and 01 arise physically from elec- 
trostatic repulsion and counterion screening, but do 



not include steric interactions due to compression of 
overlapping chains . 



VI. CONCLUSIONS 

To summarize, we have applied second-order per- 
turbation (linear response) theory to model solutions 
of spherical polyelectrolyte star-branched and micro- 
gel macroions. The theory predicts the counterion 
density profiles inside and outside of the macroions, 
effective interactions between pairs of macroions, 
and a one-body, density-dependent, volume energy 
that contributes to the total free energy of the sys- 
tem. The main conclusions are: (I) Counterions 
penetrate stars more easily than they do microgels. 
(2) Inside a star macroion, the density profile of 
mobile counterions varies more gradually than the 
macroion monomer density profile, diverging loga- 
rithmically toward the center. (3) The fraction of 
counterions trapped inside a macroion depends only 
on the ratio of macroion radius to Debye screening 
length and increases monotonically with this ratio. 
(4) Counterion screening significantly weakens the 
bare electrostatic pair interactions, which remain 
bounded up to complete overlap of macroions. (5) 
The effective pair interactions are softly repulsive for 
overlapping macroions and Yukawa in form for sepa- 
rated macroions, with amplitudes depending on the 
type of macroion. 

It is important to point out some limitations of 
the theory. First, the linear response approxima- 
tion limits applicability of the theory to dilute so- 
lutions of weakly charged macroions. The quan- 
titative range of validity depends on the relative 
magnitudes of nonlinear corrections, including three- 
body and high-order interactions, in the perturba- 
tion expansion. The same techniques that have been 
used to analyze nonlinear response in charged col- 
loids [23 can be applied to polyelectrolytes. Sec- 
ond, the mean-field approach taken here ignores fluc- 
tuations in the counterion distribution, which may 
be especially relevant for short-range interactions 
and multivalent counterions. Third, the neglect of 
chain flexibility restricts the theory to nonoverlap- 
ping macroions. This restriction may be reasonable 
for dilute solutions of sparsely separated macroions. 
However, for a sufficient concentration of macroions 
in a good solvent, chain elasticity and entropy must 
play a role. A unification of linear response theory 
and the variational theory of ref. |28l | may then prove 
fruitful. 

In principle, the predicted counterion profiles 
could be probed experimentally, e.g., by neutron 
scattering, using isotopic labelling to contrast the 
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PE chains and counterions. The macroion-macroion 
interactions may be less accessible to experiment. 
Conceivably, the solvent quality might be tuned to 
minimize the second virial coefficient between neu- 
tral monomers of overlapping macroions, effectively 
highlighting electrostatic interactions by masking 
any steric interactions. Comparisons of predicted 
and observed macroscopic properties will provide the 
most practical, if indirect, tests of the theory. Future 
applications will examine thermodynamic phase be- 
havior, especially possible implications of the volume 
energy for the stability and structure of deionized 
solutions 0. 
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APPENDIX: INTERACTIONS BETWEEN 
OVERLAPPING MACROIONS 

The bare Coulomb interaction between a pair of 
macroions at center-center separation r is given in 
general by 



Vmni{r) = — / dr' / dr' 



Pmon 



(A.l) 

For nonoverlapping macroions, spherical symmetry 
reduces the interaction to 



2^2 



r > 2a. 



(A.2) 



For overlapping macroions, the six-dimensional in- 
tegral in Eq. (|A.ip may be reduced, by exploit- 
ing cylindrical symmetry and Gauss's law, to two- 
dimensional integrals, which in turn may be evalu- 
ated analytically. For star macroions, the result may 
be expressed piece- wise as follows: 
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For computational purposes, the remaining inte- 
grals may be expressed as convergent series: 



(- 

Va 

In 



+ -ln 

a \a 



Z^e^ (9 7r 1 
2ea \2 ^ 4a ~ 2 
a — r\ r / 1 
a J ^ 2^ V2 



In 



In 



ir/ar-ir/{a~r)y 



t^mm (^) 



(r/(a- r)) 



n=l 
2n-l 



- 2 



(2n- 1) 



Z^e^ (9 7r 1 



, r < a/2 (A.5) 




'^mm (^) 



a/2 < r < a 



z^e^ rg 7r 1 r 



(A.6) 



2ea I 2 4 a 2 

' r — a 

In ( - I I In I 

a 



( a 
V r 



1 - 



a 



a 



-y 



(a/r)" -(1- a/r)' 



a < r < 2a 



(A.7) 



For microgel macroions, the bare interaction may 
be expressed more compactly as 
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Finally, the effective pair interaction between mi- 
crogels, Vcs{r), is the sum of Eq. HA.8|I and the in- 
duced interaction, obtained by Fourier transforming 



Eq. 103): 
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FIG. 1: (a) Counterion number density profiles [from 
Eqs. lHOJ, lEZJ, and lESJ] inside and outside of 

polyelectrolyte star and microgel macroions of diameter 
a = 100 nm, valence Z = 100, and effective volume frac- 
tion 77 = 0.01, in water at room temperature (As = 0.714 
nm) . The result for a hard-sphere macroion is shown for 
comparison, (b) Comparison of counterion and monomer 
density profiles for a star macroion on a log-log scale. 



12 



1.0 



O 

LL 0.5 



O 

a5 



o 
O 



star 


(a) 


microgel 

/ y 
/ y 
/ y 


^ — ■ 


/ / 
/ / 
/ / 
/ / 
/ / 
/ / 
/ / 
// 
// 
// 
A 





5 10 

Screening Constant (kg) 



1.0 



o 

CO 



LL 0.5 



O 

a5 



o 
O 



linear response (b) 
two-pliase approx. 




0.1 0.2 0.3 

Volume Fraction (r|) 



0.4 



FIG. 2: (a) Fractions of counterions [from Eqs. H25^ and 
Moi l inside polyelectrolyte star and microgel macroions 
vs. Debye screening constant k. (b) Fraction of counteri- 
ons inside a uniformly charged spherical macroion vs. ef- 
fective macroion volume fraction as predicted by linear 
response theory (solid curves) and by the two-phase ap- 
proximation of Oosawa Q (dashed curves). For each 
case, the bottom curve corresponds to ZXs/a = 8 and 
the top curve to ZXs/a = 16. 




FIG. 3: Electrostatic interactions [from Eqs. H29|l and 
I43l l1 between pairs of polyelectrolyte stars (a) and mi- 
crogels (b). Dashed curves: bare interaction. Solid 
curves: effective (bare + induced) interaction. Parame- 
ters are the same as in Fig.0 Beyond overlap (r ja > 1), 
the interaction is Yukawa in form. For overlapping 
macroions (rjcy < 1), the soft repulsion remains finite 
at complete overlap {r — 0). 
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FIG. 4: Effective electrostatic forces between pairs of 
polyelectrolyte stars (a) and microgels (b), correspond- 
ing to the interaction potentials in Fig. |^ Parameters 
are the same as in Fig. 
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FIG. 5: Macroion-size-dependent amplitude of Yukawa 
eflective electrostatic interactions [Eqs. 12911 . (I30II . and 
1431 11 between pairs of nonoverlapping stars, microgels, 
and hard spheres vs. Debye screening constant, normal- 
ized to unity at Ka = 0. 



